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The problem of infrared divergence of the effective electromagnetic field at finite temperature 
(T) is revisited. A model of single spatially localized electron interacting with thermal photons 
is considered in the limit T — > using two different regularization schemes. The first is based 
on the shift iO — > ie of the electron propagator pole in the complex energy plane, and is used to 
explicitly calculate the effective field in the one-loop approximation. We show that the matrix- 
valued imaginary part of the electron self-energy can be consistently related to the pole shift, and 
that the presence of the heat bath leads to appearance of an effective s ~ T, thus providing a natural 
infrared regulator of the theory. We find that the one-loop effective Coulomb field calculated using 
this e vanishes. The other scheme combines an infrared momentum cutoff with smearing of the 
5-functions in the interaction vertices. We prove that this regularization admits factorization of 
the infrared contributions in multi-loop diagrams, and sum the corresponding infinite series. The 
effective electromagnetic field is found to vanish in this case too. An essentially perturbative nature 
of this result is emphasized and discussed in connection with the long-range expansion of the effective 
field. 
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I. INTRODUCTION 



Interpretation of the effective (mean) field calculations is one of the most subtle questions in quantum field theory. 
It is raised every time one attempts to go beyond the realm of the scattering matrix, and is of vital importance in 
investigating such issues as the spontaneous symmetry breaking in the Electroweak theory, the quark-gluon plasma 
effects in quantum Chromodynamics, particle creation and inflation in quantum cosmology etc. In essential, this 
question consists of the following three apparently unrelated parts: I) The gauge-dependence problem, i.e., dependence 
of the effective field on the choice of the auxiliary conditions used to fix gauge invariance, 2) Measurability of the fields, 
i.e., consistency of the formal field-theoretic predictions with the principal limitations of quantum measurements, and 
3) The presence of infrared divergences. 

The first two of these particular problems are formidable in Yang-Mills theories, which is why the question of 
interpretation of the effective field is still far from its resolution there. In quantum electrodynamics, on the contrary, 
they are relatively simple and amenable to general analysis. In fact, the proof of gauge independence of the effective 
quantities such as Heisenberg-Euler Lagrangian, the effective (mean) field, etc. is not difficult and has been given 
in quite general form Furthermore, measurability of the electromagnetic field was proved in the fundamental 
works by Bohr and Rosenfeld 0, Q - In particular, it was argued in [l| that the measuring apparatus can in principle 
be designed so as to take into account radiative corrections despite the presence of ultraviolet divergences. This 
consideration was subsequently generalized and carried over to the gravitational interaction in 0] . 

What remains unsolved, however, is the problem of infrared divergences of the fields produced by elementary 
particles. Namely, it is known very well that the vertex formfactors of massive field quanta, used in constructing 
effective fields, are infrared-divergent. In this respect, electromagnetic interaction has many similarities with gravity, 
in particular, regarding resolution of the infrared catastrophe. Unfortunately, the corresponding proof [a, 0] of 
cancelation of the infrared divergences in the S-matrix does not apply to the effective fields. In fact, the very statement 
of the problem in the two cases is quite different. The Bloch-Nordsieck theorem states that infrared divergences in 
the radiative corrections to the scattering cross sections exactly cancel those due to the emission of real soft photons. 
At the same time, the latter do not appear in the effective field formalism at all. An opinion can be found in the 
literature (see, e^o., 0) that this problem can be overcome by using the effective potential defined via the S-matrix 
according to [H, |sj]. It is known, however, that such definition is highly ambiguous 10]: There are infinitely many 
potentials which lead to the same S-matrix, and not all of this arbitrariness can be traced to the gauge freedom in the 
theory. As a result, even to achieve agreement with the classical theory requires introduction of certain prescriptions 
at each order of the perturbation expansion. Apart from this difficulty, such a definition would not help investigating 
the issues mentioned at the beginning. 

By this reason, the effective field formalism can be applied, strictly speaking, only to the fields produced by classical 
sources. If, for instance, the mass of a charged particle is sufficiently large, radiative corrections to its interaction 
with the electromagnetic field can be neglected, thus putting the question about their divergence aside. In fact, 
assumptions of this kind underlie the classic calculation [TTj and its generalizations [T^ |. However, this limitation 
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looks rather unsatisfactory from the theoretical point of view, especially in the light of the positively resolved issue 
2) above. 

This problem is sharpened at finite temperatures, because the heat-bath effects introduce new low-energy singu- 
larities to the photon propagator. These additional singularities are generally believed to worsen infrared properties 
of the Feynman integrals. Previous investigations of the infrared problem have been aimed mainly at generalizing 
the Bloch-Nordsieck and Kinoshita-Lee-Nauenberg theorems to nonzero temperatures [3, EH EH, EH, E3, EH, EH]- 
It should be mentioned in this connection that extending the S-matrix formalism to finite temperatures is itself a 
nontrivial task. Some general results in this direction were obtained in [16[, where an S-matrix connecting states with 
different temperature was constructed. An independent definition of the finite-temperature S-matrix is given in [ItJ 
where infrared finitencss of the S-matrix is proved in the eikonal approximation. These results were subsequently 
improved in Though in a different way, cancelation of the lowest-order infrared divergences at finite tempera- 
ture was earlier demonstrated in (l3j . However, it was argued later that the infrared divergences treated either way 
still persist at arbitrary temperatures It is important to stress that all these developments proceed within the 
conventional approach based on summation of the scattering cross sections over suitable classes of processes. As was 
already mentioned, their results do not apply to the effective fields. 

There is an argument often used in the literature, according to which infrared divergences in the formfactors can 
be discarded within the long-range expansion of the electromagnetic or gravitational fields. Indeed, it is known 
that the infrared-divergent contributions considered as functions of the momentum transfer, p (which is the variable 
Fourier-conjugated to the distance from the field source), are analytic at p — 0. Though infinite, such contributions 
consequently give rise to terms proportional to the Dirac delta-function or its derivatives in the coordinate space. 
The divergent part of the effective field thus turns out to be confined to the region of source localization. Yet, this 
argument cannot be recognized convincing. The point is that such terms appear at every order of the perturbation 
theory, and each order adds new derivatives acting on the delta-function. Therefore, the question as to locality of the 
divergent contribution can be settled only after summing the infinite series of differential operators, because the sum 
of local operators may well turn out to be a nonlocal operator. 

Turning back to the infrared problem in the general case, we would like to stress that inclusion of the finite 
temperature effects into consideration is a question of principle: no matter how small temperature might be, it is 
never zero exactly, and the question as to whether T can be replaced by T — can be answered only by 
investigating the general nonzero-temperature case. In other words, continuity of this sort, being a necessary physical 
requirement for the very possibility to neglect the heat-bath effects, is to be proved rather than postulated. It may 
well turn out that in view of the aforementioned principal inequality T ^ 0, the notion of the vacuum effective field 
is to be defined as the limit of its finite-temperature counterpart for T — > 0. The purpose of the present paper is to 
examine the infrared problem from this point of view. In this respect, our investigation differs from the standard 
approach based on the high-temperature expansion, which is commonly used in actual calculations (consideration of 
the low-temperature case can be found, e.g., in [lj, EH])- We restrict ourselves to the electromagnetic interaction 
because its diagrammatics is much simpler than gravitational, and admits in certain cases partial summation of the 
perturbation series, which cannot be performed in the presence of self-interaction of the gauge field. We introduce a 
simplified model for a single spatially-localized charge interacting with the heat-bath photons at equilibrium, which 
neglects the influence of the photon interactions on their statistical distribution. This simplification is in fact a valid 
approximation in investigating the infrared problem, because it is the low-energy asymptotic of the photon distribution 
that is only important in this case, while the radiative corrections to this distribution vanish at zero momentum. This 
model requires some modification of the diagrammatic rules, which is described in Sec. [TTJ To examine the infrared 
divergences of the effective electromagnetic field, we use two different regularization schemes. The first, called e- 
regularization for brevity, consists in the usual finite shift iO — > ie of the poles of the electron propagator in the 
complex energy plane. This regularization is quite natural in view of the fact that the temperature effects lead 
to appearance of an imaginary part in the electron self-energy, thus providing an effective infrared regulator which 
replaces e by a quantity ~ T. This scheme is also very convenient in actual calculations as it reveals some important 
features hidden in other regularizations, such as the cancelation of counterterm contributions to the effective field, 
considered in Sec. IIII1 However, it does not admit factorization of many-loop infrared contributions, which hinders 
partial summation of the perturbation series. The other scheme, called A-regularization in what follows, is based 
on the introduction of an infrared cutoff, and is described in Sec. IIV Al This scheme does admit summation of the 
infrared contributions similar to that employed in the standard treatment of the infrared catastrophe. This summation 
is performed in Sec. IIV Bl The general conclusion drawn in Sec. [V] is that both regularizations suggest vanishing of the 
effective electromagnetic field calculated using the perturbation theory, thus revealing an essentially non-perturbative 
nature of the infrared problem. The paper has an Appendix where some technical aspects of our work are discussed 
on the basis of the Ward identities. 
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II. PRELIMINARIES 



Consider the electromagnetic field produced by an electron of mass m, which is on average at rest and interacts 
with the virtual as well as real photons in equilibrium at finite temperature 1 T <C m. Neglecting the effects related 
to the quantum spreading of the electron wavefunction, the field is purely electric, so that it is sufficient to find the 
effective scalar potential 



Af(x) = H^Tr (A {x)e- pH *>g) , N = Tr (e- pH *g) , P = l/T. 



(1) 



Here Aq(x) is the Heisenberg picture operator of the scalar potential, the Hamiltonian of free photons, g the 
electron density matrix, and the trace is over all photon states as well as the single electron states. The matrix g is 
assumed to describe a non-thermalized state, viz., a state in which the electron is spatially localized. Having chosen 
Hcf, to be the free photon Hamiltonian we thereby omit corrections to the photon distribution due to their interactions. 
It is readily seen that within the one-loop approximation, these interactions do not change the effective field indeed. In 
fact, in investigating the infrared divergence problem at higher orders they can be neglected as well. This is because in 
this case, it is the low-energy asymptotic of the photon distribution that only matters, while the radiative corrections 
to the photon Green functions vanish at zero momentum to all orders. Also, the electron Hamiltonian might have 
been added to the exponent in Eq. (TTJ). However, since there is only one nonrelativistic electron, and T <C m, its 
contribution can be easily shown to factorize and cancel the corresponding contribution to the normalization factor 
N, leaving the value of the effective field unchanged. 

To evaluate the right hand side of Eq. |T]) , we shall use the general framework of the real time approach [20| ■ Since 
there is only one charged particle present which is not in thermal equilibrium with the photon bath, the standard 
finite-temperature-field-theory techniques are not directly applicable to the calculation of Aq R . The relevant extension 
of the real time formulation can be obtained as follows. According to the Schwinger-Keldysh method [2l|,l22|], the 
expectation value (|TJ) can be written as 



exp i 




A (x) 



(2) 



where Lj is the interaction Lagrangian, and the ^-integration is along the contour C — C\ U C2 in the complex time 
plane, shown in Fig.[TJa). The subscript in indicates that the operators are taken in the interaction picture; T c is the 
usual ordering along the contour C. Concerning the choice of the time contour, the following circumstance must be 
mentioned from the outset. It is a common procedure of the real time techniques that the factor e~@ H , with H the full 
Hamiltonian, is promoted into an integral along the interval [ti, ti—i0\ in the complex time plane, followed by a contour 
deformation shown in Fig. [T^b). Taking the limit ti — > — 00, tf — > 00 then allows one to get rid of the contributions 
of the contour segments C3 or C4. This procedure presupposes eventual factorization of such contributions, resulting 
from the vanishing of the two-point Green functions having one argument on C3 or C4 and the other on C\ or C2. It 
is known, however, that in theories containing massless particles, this transformation is jeopardized by the infrared 
singularities [2(j. It turns out that the model we consider is the case where the trick with the contour deformation 
is not legitimate. Specifically, an explicit calculation shows that the one-loop diagrams involving the Green functions 
which connect with C3 actually do not vanish. In particular, omitting them makes the effective electromagnetic 
field complex- valued. On the other hand, the replacement of H by the free photon Hamiltonian makes it unnecessary 
to follow the above-mentioned procedure. Thus, to avoid the complications associated with the contour deformations, 
we adhere to the original Schwinger-Keldysh formulation. 

Perturbation expansion of the right hand side of Eq. ([2]) leads to the sum of terms of the form 



(3) 



A J d*j/...T*{e-fi H *T c [Mz)Mv)MM) 
xTr (T c [: ^(y)rHy) ■■ i>{v')l l W) : Ao) ■ 



Hereon all the field operators are taken in the interaction picture, and the subscript in is omitted. As usual, the 
particle propagators acquire 2x2 matrix structure as the result of T c -ordering. By virtue of the Wick theorem, the 
first factor in the above expression is represented as the sum of products of various finite-temperature photon two- 
point functions, while the second a similar sum constructed from the vacuum Green functions for the electron-positron 



1 We use relativistic units h = c = 1. Also, the Minkowski metric is ijfu, = diag{+l, — 1, — 1, — 1}. 
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FIG. 1: (a) Integration contour used in Eq. ((2}. (b) The contours commonly used in the real-time techniques, parameterized 
by a 6 (0, 1). The limit ti — > —oo,tj — > oo is understood. 



field, and the sing le factor Tr (ip(y')ip(y")g). Each propagator bears indices 1 or 2 assigned to its ends, according to 
whether the end belongs to the contour branch C\ or C2, and connects vertices of the corresponding type, with an 
extra minus sign for each 2-type vertex. The observation point x is assigned index 1. In momentum space, the matrix 
photon propagator in the Feynman gauge reads D^{k) = rj^D^\k), where 



D^\k) = -b^ 22 \k) = 



- 2mn{k)8{k 2 ) 



k 2 +i0 

D [21) (k) = -2iri[e(k ) + n{k)}6(k 2 ), 
D (12 \k) = -2iri[9(-k Q )+n(k)}5(k 2 ), n(k) 

while the electron matrix propagator D e {k) has the elements 



D™{k) = -D^(k) = / + W - n , 

m — fc — iv 

D (21) (k) = 2ni8(k 2 -rn 2 )6»(fc )(^ + m), 
D^ 2 \k) = -2iriS(k 2 - m 2 )6»(-fc )(^ + m) , 



(4) 

(5) 
(G) 



where the tilde symbolizes the special operation of complex conjugation which treats the 7-matrices as real quantities: 
7 = 7. It will be convenient to represent S e (fc) in the form 



D e (fc) = M(k) 

M(k) 

which is established directly using Eqs. dU)-©. 



Dp 
-Dp 

( 1 -e(-ko) 
V 0(ko) 1 



M(k), D F {k) 



1 



k 2 - i0' 



(7) 



III. ONE-LOOP CONTRIBUTION TO THE EFFECTIVE FIELD 



A. e-regularization 

Infrared singularities in the electromagnetic formfactors and self-energy contributions to the effective field, as well 
as ultraviolet divergences require introduction of intermediate regularization. The simplest and very convenient choice 
we make in order to investigate the effective field in the one-loop approximation is the usual finite shift of the electron 
propagator poles in the complex A:°-plane. Namely, we replace the function Dp(k) in Eq. ([7]) by 

D F (k) = 2 * e>0 

m A — k z — le 
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(similar regularization of the photon propagator is unnecessary) . Accordingly, the 5-function appearing in the electron 
propagator is regularized as 

2m8(k 2 - m 2 ) - ' ' 



m 2 — k 2 — is m 2 — k 2 + is 

Obviously, e plays the role of infrared regulator, and for brevity, we call this scheme e-regularization in what follows. 
In fact, the introduction of e-regulator is necessary in the real-time formulation of thermal field theory in order to 
correctly treat the products of singular functions appearing in the Feynman diagrams. As we shall see, it helps revealing 
some nontrivial features of the one- loop contribution to the effective field, which are hidden in other regularizations. 
To make the e-regularization physically meaningful, one should impose some restriction on the value of e with respect 
to the physical parameters of the theory. Since e has dimension [mass] 2 , such a restriction should be smallness of e 
in comparison with some quadratic combination of T and m. We assume throughout that e satisfies 

e < vaT. (8) 

It will be seen below that this condition appears naturally in the theory. In addition to that, we introduce a momentum 
threshold, A, separating the infrared effects. It satisfies T < A < m, and identifies the photons with k < A as "soft." 
As to the usual ultraviolet divergences, they are supposed to be regularized using some conventional means, say, the 
dimensional technique. 

Upon transition to the momentum space, Aq S (x) takes the form 

A f{x)=Y,Jj ^-^e lpX D {11) (p)Q^(q,q+p)^(q + p)J(p,q)u a (q), (9) 

where the matrix-valued function J(p, q) is essentially a sum of products of several momentum-space propaga- 
tors. In this formula, the values of q° and p° are fixed by the mass-shell condition q° — q 2 + m 2 > and the 
energy-momentum conservation. Since the field-producing charge is non-relativistic, p° = (q + p) 2 /2m — q 2 /2m . 
The bispinor amplitudes and the momentum-space density matrix are normalized on unity: 

u a u a =1, J J^j3 Q<t<t(Q> l) = 1 • 

It is easily verified that in the tree approximation, 

J(p, g ) = -e 7 . (10) 
The one-loop diagrams representing Aq S are shown in Figs. HHU 




FIG. 2: Vertex contribution to the effective field. 



It is not difficult to see that the diagrams (c) and (d) in Fig. [2] are irrelevant. This is because they involve the 

(12) 

function D e (q — k) ~ 9(ko — qo), and therefore do not contribute at small loop momenta. By this reason, diagrams 
involving this function are omitted in Figs. [3J [4j Less trivial is the fact that the counterterm diagrams do not contribute 
either (the procedure of introducing temperature-dependent counterterms is discussed in detail e.g., in (25l. l26l|). For 
the proof, we note that the vertices in Fig. [H generated by the counterterm Lagrangian, are contracted with the 
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FIG. 3: Electron self-energy contribution to the effective field. 




FIG. 4: The self-energy counterterm diagrams. 



bispinor amplitudes on the mass shell. Therefore, the corresponding vertex factor is proportional to the unit matrix. 
Denoting l c the proportionality coefficient, we find for the sum of diagrams in Fig. [4] 



m) —-- ^ + i 



IE \ —IE IE I —IE 



= 0, (11) 



since (^+^)7° — > 7717°, 7°^ — > 7717° upon contraction with u,u on the mass shell. We can interpret this fact by saying 
that when calculating the effective field, the usual renormalization prescriptions for the electron propagator can be 
assumed fulfilled automatically. Thus, only diagrams in Fig. (Ha)lUb) and Fig. [3Ja)-[31[d) remain to be considered. 

Our primary aim below is to extract the small loop- momentum parts of the above diagrams, and to show that the 
one-loop contribution to the effective field is infrared- finite. In doing this, we consider arbitrarily small but finite 
temperatures, and calculate asymptotics of the diagrams as T — > 0. The calculation reveals two types of infrared- 
divergent terms in these asymptotics. Their leading at T — > terms are proportional to ln(/3e), while the subleading 
ones are linear in temperature. It turns out that the singular part lne of the former exactly cancels the corresponding 
vacuum contribution, so that the finite remainder proportional to In f3 determines the low-temperature behavior of 
the vertex and the electron self-energy contributions to the effective field. The role of the other terms is different. 
The point is that the linear in T contribution to the electron self-energy contains an imaginary part, and acts as a 
natural infrared regulator of the theory. In other words, taking into account the subleading terms is equivalent to 
replacing e — > const • T. It turns out that this gives rise to some finite temperature-independent terms in the vertex 
formfactor which exactly cancel the tree contribution. 



B. Electron self-energy contribution 



To begin with, let us show that the one-loop contribution to the normalization factor Ji vanishes. This is conveniently 
done using the following representation for the matrix electron self-energy [23| 



=W = «-'(«> I 0' -S, ) - ( - %o)( E F - IS,) ' -X I ' < 12 » 
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FIG. 5: One-loop contribution to the normalization factor 3\f. 



where M is dchncd in Eq. ([7]), and Sf = S'* 11 -' (g)- It must be mentioned that this representation is generally valid 

fll) f22) (12) (21) 

only in the limit e — > 0. This is because its derivation is based on the identity T) e + D e = T> e + D e for the 
exact electron propagator, which is a direct consequence of the time-ordering prescription, and hence, is affected by 
the e-regularization. However, explicit calculations show that it holds true also for the leading terms we are interested 
in, and so we use Eq. (|12p below to render the formulas more compact. The one-loop contribution to N, shown in 
Fig. [5j reads 



d 3 q 



<T,<T f 



Using Eq. p2[) shows immediately that the integrand is zero identically (in this case, the e-regularization is unneces- 
sary) . 

Moreover, it turns out that the real Lorentz-invariant parts of the electron self-energy (reality being understood in 
the sense of the operation ~ ) also cancel upon substitution to the effective field. To see this, let us write down the 
contributions of diagrams in Fig. [3] to the function J(p, q) 



1 



J(3a) = ie'E F (q + p)-($+ft + m)j° 



'(3b) 



'(3c) 



J, 



(3d) 



ie {£ F (q + p) - Z F {q + p) } J (jf + f + m) 7 ° 
iej°(d + to)-Ef(<?) 

r £ 



(13) 
(14) 
(15) 
(16) 



Recalling that these expressions are sandwiched between bispinor amplitudes, we may replace <b + f — > m in the 
first three of these expressions, and <b — ► m in the last one. Taking into account also the mass-shell conditions 
q 2 = (q + p) = m 2 , we conclude that the Lorentz-invariant part of their sum involves only the difference Sjr — Hp. 
It follows that the vacuum contribution is 0(e), since it is Lorentz-invariant and real for e = 0. 
The explicit expression for Yip at the one-loop level reads 



d A k 7^ + $ + m) 7/1 
(2tt) 4 m 2 -(q + k) 2 - ie 



1 



k 2 + i0 



2mn{k)8{k 2 ) 



(17) 



In view of what was just said, the ultraviolet divergence of Si? is of no importance for the present computation. 
Furthermore, since its imaginary part comes from integration over small k, we may cut the integral at k = A 3> T. As 
T is arbitrarily small, A can be taken small enough to neglect $ in the numerator as well as k 2 in the denominator. 
Then the vacuum contribution (corresponding to the first term in the square brackets) takes the form [we have to 
retain the 0(e:)-terms because they are combined with the mass-shell electron propagators in Fig. [31 which are 0(l/e)] 



Svac 
F 



[a) 



-ie 2 (2m — d) 



2tt 2 



(2m - , 



d^k 



1 



(2tt) 4 2kq + ie k 2 

A 

dQ f ndn 
4ir / 2na — ie 



i0 



-ig^(2m-. 



In- 



where a = kq, k = \k\ = 1, J dil is the integration over directions of fc, and "• • ■ " stands for real terms, or terms of 
higher order in T, e. 

To find the heat-bath contribution, let us use the condition ([8]) to introduce an auxiliary intermediate scale A such 
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that e/m < A < T. Then the loop integral can be evaluated as 

oo 



1 



2tt 2 

^2 



(2m - . 



47T 



— \n!— + — / 

m m p J 4-7T a 



1 k\ / 1 1 

/3 2 J \2na — ie 2na + is 

+ ■•• , 



2a 2 ./ k 

A 



Adding the two contributions yields 



4ln(/3m) + ^(a- 1 } 



where the angular brackets denote averaging over directions of the 3-momentum, 
Using the condition |g| < m, this expression can be further simplified to 



6 



— ln(pm) + — 
m z pm 



(18) 



which will be sufficient for our purposes below. We note, first of all, that the heat-bath contribution proportional 
to lne cancels that of the vacuum part. This fact allows us to use Hp to construct a function of e, regular in a 
vicinity of e — 0. Indeed, T,p was defined initially for e > 0, and the found cancelation of logarithmic singularities 
means that it can be continued analytically to the region s ^ 0, and that this continuation is unique. Let us denote 
it £^ = T<p(q, s). We need this continuation because Eq. (fTgj) shows that at finite temperature, the electron mass 
acquires an imaginary part ~ aT (a — e 2 /4ir is the fine structure constant), and this part is positive. 2 This means in 
turn that its mass squared becomes m 2 + iE, where m is the (renormalized) vacuum mass, and E ~ amT. In other 
words, the finite-temperature effects introduce natural infrared regulator E, replacing in effect e — > e — E, so that 
removal of the infrared regularization is now to be understood as e — > — E in the analytically continued functions. 
Specifically, the electron self-energy becomes in this limit (g, — E), while the mass-shell electron propagator, 
(f + m)/(iE). 

In order to determine the precise value of E, we have to relate the matrix-valued imaginary part of the electron 
self-energy to the scalar shift of the propagator pole. Since a change of E changes the electron self-energy, which 
in turn affects E itself, one has to use Eq. (fl8|) to find the fixed point for E, i.e., the value which upon replacing 
e — > — E in the electron propagator would give the same E as read off from the electron self-energy calculated using 
this propagator. Let us write, for brevity, ImEf(g) = ioxfy + ia2m, where ai t 2 — 0(a) are real quantities whose 
explicit expressions can be found by replacing e — > — E in Eq. (TJ5J) . Using the fact that ReT,p(q) vanishes on the 
mass shell (this condition is satisfied automatically, Cf. Sec. IIII A|) . we have for the electron propagator near its pole 



1 



+ m + ima-i — Ho\ 



q 2 + 2i(q 2 o\ + m 2 (T2) + q 2 o 2 — m 2 a 2 



Since the leading terms in <j\ t 2 are O(a), and do not contain ln(/3m), within the accuracy of our calculation this 
expression reduces on the mass shell (q 2 = m 2 ) to 

$ + m 



2im 2 io\ + (72 ) 



2 The fact that the mass eigenvalues in thermal field theory can have positive as well as negative imaginary part was emphasized in 

am. 
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On the other hand, according to the definition of E, the mass-shell electron propagator is 

<j + m 



lE 



An 2 



—n ln(/3m) + — 
m z pm 



Comparison of the two formulas yields E — 2m 2 (<7i +(73), or, using Eq. (jTH 

E = 

Thus, we find 

E = 



2ll(3 



1 + ^HPm) 



(19) 



It should be emphasized that the temperature in this calculation does not have to satisfy e 2 ln(/3m) < 1, so that 
Eq. ([19]) is valid for all T down to zero. Finally, using this result in the sum of the four contributions ([13]) - (|T6|) . we 
arrive at the following expression for the electron self-energy contribution to the function J(p, q) 



Because of the mixing of orders caused by E = 0(e 2 ) in the denominator, this contribution turns out to be a first-order 
quantity. 



C. Vertex contribution and nullification of the one-loop Coulomb potential 

At last, let us consider diagrams in Fig. [5] To extract the leading term as T — > 0, we again introduce the cutoff A, 
and replace m 2 — (q — k) 2 — ► 2kq in the mass-shell electron propagators. Performing the standard 7-matrix algebra 
with the help of the formula (d + m) r ) tl u cr (q) = 2q fl u (7 (q), one then finds that the infrared contribution of diagrams 
(a), (b) in Fig. [5] is contained in the expression 



J m = Aie 6 ^(m' -p z j2) 



(27T) 



1 



1 



k 2 + iO 2kq + is 2k(q + p) + is 



1 2mn{k)5{k 2 ) 2m8{k 2 )9{k a ) 



2kq — ie 2k{q + p) + is 2kq — is \2k(q + p) — is 2k(q + p) + is 



1 



Integration over fco yields 



e 3 7 



1\ 

(2m 2 -p 2 ) J^-JdKK 



1 



1 



2n(n) 



1 



2na — is 2nb + is 2na — is 2nb + is 



(20) 



where b = k(q +p). It is seen from this expression that in the formal limit T — > 0, the second term in the integrand 
vanishes (because n(n) — » 0), leaving us with the vacuum contribution given by the first term. However, this reasoning 
is not actually legitimate since it is drawn from consideration of the integrals which diverge for e — > 0. Taking the 
limit T — > keeping s fixed is indeed only formal as it violates condition (jSJ). At the same time, it was shown in the 
preceding section that the finite-temperature effects introduce natural scale for e, namely, E ~ amT, and this value 
automatically satisfies ©, because a < 1. Under this condition, the k- integral can be written, using the auxiliary 
intermediate scale A, 



A 

/ 



dn k 



1 



1 



2kq — is 2nb + is 



2 \ 1 1 

(3k J 2kcl — is 2Kb + is 



1 

Aab 



dn 

K 



2dK 



We see that the logarithmic singularity of the vacuum contribution is exactly canceled by the temperature-independent 
term coming from the expansion u(k) = 1/(/3k) — 1/2 + O(fiK). The other integrals are evaluated with the help of the 
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formulas 

A 

r , 1 1 ?r l l 

/ w t-* — — 

J 2na — is 2nb + is 2e a + b AXab ' 
o 

oo 

^^T r -^-^)] + - x+ H^ A^O (21) 

A 

(7 is the Euler constant). It remains to set e = —E, and use the formula (119[) . To this end, it is to be noted that 
although no restriction on the momentum transfer, p, is imposed a priori, the fact that Eq. (|19p was derived under 
the condition \q\ <C m implies that p also must satisfy \p\ <C to. Then 6 = a = to, and the leading terms in the vertex 
contribution to J(p, q) take the form 

■/(a) (P, 8) = ^(2m 2 -p 2 )HPm){(ab)- x ) 

- e7 °m (2 - 4> ) ((a + 6)- 1 ) (l + ^ M/m) = • ( 22 ) 

Using the obtained results, it is easy to demonstrate nullification of the Coulomb field of a localized charge. The 
condition \p\ <C m means that the electric field is measured far enough from the charge. More precisely, if the electron 
is localized near xq, then the distance between this point and the point of observation, r = \x — x \ 3> 1/to. Adding the 
vertex and self-energy contributions gives Jn) + J(s) — e^°. It follows that the one-loop contribution to the function 
J(p, q) exactly cancels the tree value ([TU]) . so that 

Af (x) =0 for r > 1/to. (23) 

Thus, independently of the particular choice of the electron density matrix, the leading long-range term of its electric 
potential - the Coulomb field - vanishes upon account of the one- loop radiative contributions. 

This analysis can be extended to higher orders. For instance, it is not difficult to see that two-loop diagrams 
contain terms proportional to (T/e) 2 , and therefore, contribute 0(e)-terms to the effective potential. In particular, 
they modify the quantity E, which in turn changes the structure of the 0(e 3 )-terms found at the one-loop level, etc. 
Unfortunately, e-regularization does not allow extension of the above result to all orders, because the mixing of orders 
in the coupling constant obscures the structure of the perturbation series. 



IV. HIGHER-LOOP CONTRIBUTIONS TO THE EFFECTIVE FIELD 



A. A-regularization 

In addition to mixing the perturbation series, the e-regularization has the drawback that it breaks the usual 
factorization of many- loop infrared contributions (see Sec. IIV 51 below) . To investigate their structure, we shall now 
introduce a regularization that does not lead to mixing of various orders of the perturbative expansion, and allows 
direct summation of the divergent contributions. 

First of all, we introduce the usual infrared cutoff Ao <C T, so that all loop momenta are now restricted to the interval 
Ao < k < A. Still, this cutoff is insufficient for the purpose of regularizing diagrams with the self-energy insertions. In 
order to make these and the corresponding counterterm diagrams meaningful, it is necessary to regularize the mass- 
shell propagators (see Figs. [31 H]). Leaving the mass shell does not help in this respect, as it precludes factorization of 
the infrared contributions. By this reason, we introduce the following smearing of the Dirac delta- functions expressing 
conservation of the 4-momentum in the interaction vertices 

5 4 ( qi +k-q 2 )^A x ( qi + k~q 2 ), (24) 

where A\(w) satisfies 

j d 4 wA x (w) = 1, A A H = A x (-w), (25) 

A A (w^0)^0 as A^O. (26) 
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A convenient choice is 

A,(-)-^-p(-^). m 

It will be assumed in what follows that the parameter A characterizing the width of the smeared delta- function satisfies 
A <C Ao. We shall also use the following notation 

aJM-h/M^M 2 ^)- (28) 



a) 




b) 




FIG. 6: Self-energy insertion into external line of a general diagram (a), and the corresponding counterterm diagram (b). 



Next, to ensure convergence of the effective field in the limit A — > 0, we have to introduce counterterms into 
Lagrangian. To be specific, let us consider diagrams of the type shown in Fig. [6^a). The external line with the 
self-energy insertion contributes a factor 



d A w\ / d i W2&\{w\)£i.\(w2) 



2q{w\ + w 2 ) + i0 



+ v/i+v/ 2 + m)S (11) (g + wi) 



1 



d i wA* x (w) 2qw + M (<f + yt + m) [^(g) + 0(A) 



(29) 



Here the quantity £( n )(g) is taken on the mass shell. The renormalization prescription that the propagator poles be 
at the physical mass requires vanishing of this quantity, which can be met by introducing a counterterm. The above 
expression shows that the delta-functions appearing in the two-point counterterm vertices are to be regularized as 



S A (qi - qi) -> A^(<7i - q 2 ). 



Then the counterterm diagram shown in Fig. [6](b) reads 



d 4 wA* x (w) 



1 



2qw + iO 



(^ + l/} + m)^ n \q). 



(30) 



(31) 



The remaining contribution coming from integration of the 0(A)-term in Eq. (|29[) will be taken care of below when 
summing all infrared many-loop contributions. It is worth noting that in the limit p — > 0, the counterterm contributions 
to the effective field cancel each other, just as in the case of e-regularization. This can be verified directly repeating 
the calculation of Sec. IIIII [Cf. Eq. (fTTj) ]. However, this property is violated at finite p by the Lorentz-noninvariant 
integration in Eq. (f3"Tj) . 



B. Factorization of infrared contributions 



Now that the finite-momenta contributions to the electron self-energy have been canceled by counterterms as 
discussed in the preceding section, it remains to take into account contributions of small virtual photon momenta. 
The general structure of diagrams to be considered is shown on Fig. [7J In terms of the function J{p, q), the sum of 
all such diagrams can be written as a series 

oc 

J(p, q) = -e7 / A (p, q) £ (e 2 ) N I N (p, q, A), (32) 

JV=0 
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where the factor I\(p,q) corresponds to photons with momenta k > A. Let us first show that the £-regularization 
does not admit the usual factorization of infrared contributions. In has the following general form 



N N-m 



I N {p,q,A) = w 



d k\ d kyiQ 
(2^)4 "' (2^ 



m 2 (™*+"/)( m 2 -p 2 /2) n °D( u \h) ■ ■ ■ D {11 \k no ) 



m—0 Uf—O 

F^(q, k u ..., k no )F^(q+p, k u ..., k„ 
nj!n/!n ! 



(33) 



where N is the number of virtual photon lines, of which (rif) reside on the ingoing (outgoing) electron line, while 
the remaining n = N — rii — rif connect the two electron lines. In the case of e-regularization, the function F™ can 
be conveniently written as 



F„"° (q, k u . . . , k no ) =J d V • • • d A w 2n+na D (1 V (w 2 )---D^ (w 2n ) 



xS(k! - W 2n +l) ■ ■ -S(k no - W 2n+no )S(w 2 + Wl) • ■ ■ S(w 2n + W2n-l) 
1 1 1 



E 



wiq + ie (wi + w 2 )q + is (w\ H h w 2n+no )q + is 



+2-terms, 



(34) 



where the sum is over all permutations of indices 1, 2, . . . , 2n + no, and "2-terms" denotes the contribution of diagrams 
involving 2- vertices. A partial permutation with respect to the indices 1,2 of the expression in the square brackets 
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yields 



1 



1 



wiq + ie w 2 q + ie 



(w 1 +w 2 )q + iej (w% H lw 2n+no )q + ie 



(35) 



It is seen that the usual factorization is hindered by the term proportional to e, which is not negligible because 
(wi + w 2 )q may take on values as small as e. Things are different in the A-regularization which replaces Eq. (|34jl by 
the following 



d 4 k no +i ■ ■ ■ d A k 



n +n 



d 4 Wl ■ ■ ■ d 4 w 2n+no D^ (fc no+1 ) ■ • • (k no+n ) 



xA(w 2n+1 - ki) ■ ■ ■ A(w 2n+no - k no ) 

xA(wi - k na+ i) ■ ■ ■ A(w 2n -i - k na+n )A(w 2 + k na+1 ) ■ ■ ■ A(w 2n + k no+n ) 

1 1 1 
2 — - ; — — : : : ; — — — + 2-terms. 



w\q + ie (wi + w 2 )q + ie (tui H h w 2 „+„ )g + ie 



(36) 



Now the second term in the large parentheses in Eq. ([33)1 is O(s/Xo), and it is not difficult to check that so are all 
other terms proportional to e, appearing in the course of the permutation. Therefore, under the assumption mAo ^> e, 
they can be omitted, and we may consistently set e — 0. 

Next, let us show that the 2-terms do not contribute to the effective field in the A-regularization. Figure [8] depicts 
general structure of the corresponding diagrams. By the reason explained at the beginning of Sec. IIII1 diagrams with 
a 1-vertex appearing to the left of a 2-vertex vanish. If a diagram has m > 1 vertices of the 2-type on the outgoing 
electron line, then the corresponding sum over permutations reads 



E 

perm 



1 



1 



w±q - iO (wi + • • • + w m -i)<? — iO \ (wi + 



1 



iO 



1 



1 



(wi + • ■ ■ + w m )q — iO ) (wi + ■ • • + w rn+ i)q + iO (wi + 
1111 



y- 



perm 

0, 



wiq — iO w m q — iO Wiq — iO w m q — iO J (w% + ■ ■ ■ + w rn +i)q + iO 



W2n f +n )q + iO 
1 



where we explicitly performed the permutation over indices 1, m, and used the factorization rule proved above. In 
the case m = 1, expression in the square brackets takes the form 



1 



1 



1 



1 



w±q + iO wiq —iOJ (wi + w 2 )q + iO (wi + ■ • • + w 2nj+no )q + iO 

This vanishes, too, because Wi ^ in the A-regularization. Thus, the contribution of the 2-terms is zero. This result 
is discussed from a more general point of view in the Appendix. 

We are ready to perform summation of the infrared contributions in the series (|32p . This is very similar to the 
standard treatment of the loop infrared divergences (see, e.g., [13|). The sum over permutations in Eq. (|36p becomes 



1 



1 



1 



wxq + iO w 2 q + iO w 2n+na q + iO 
Substituting Eq. ([36]) into Eq. (l33|) . and designating 



(37) 



-iq s qt 



d 4 k 

(27T)* 



d wi id w 2 



A(wi - r) s k)A(w 2 + rj t k)D^(k) 



(wiq s + i0)(w 2 q t + iO) 



where s, t = 1, 2, qi = q, q 2 = q + p, rji = 1, r/ 2 = — 1, we obtain 



I N (p,q,A) =J2J2 



d ll d 22 d 12 _ ( d 



d 22 + 2d 12 ) N 



m=0 n f =Q 1 ' u 



^12^ 



(38) 



Preliminary version - October 30, 2009 



14 



Under the assumption A <C Ao, d s t can be evaluated as 

d 4 k 



(27T) 



d st = iVsVtq s qt 
Therefore, in the limit of removed A-regularization, 

d st = iVsiltqsqt 
Substitution of Eq. (38|) into Eq. ((32j) gives 



<W^ 2 (« 



o 



d 4 k D^jk) 
(2tt) 4 (kq s )(kq t ) 



J(p, q) 



-cy exp 



-(d 11 +d 22 + 2d 12 ) 4(p,9) 



The integrals d s t can be computed in the same way as the J(2)-integral in Sec. MI Cl One has 



, VsVt 



f d&l 1 f dn r „ 

/ ; / — l + 

y 4^ (fc 9s )(fcg t ) y « 



and the use of Eq. (|21[) brings this function to the form 



, VsVt 
st = ~4^ qaqt 



dn 



4tt (kq s )(kq t ) 



/3A 



+ ln(/3A) + (7 - ln(27r)) 



Thus, retaining the singular terms in the exponent gives finally 

J(p,q) = -e 7 °exp (- [l - (m 2 -p 2 /2) ((ab)- 1 )} 



pA 



^"a(p, q) ■ 



The expression in the exponent is negative for all p. For small momentum transfer, for instance, one has 



J(P,Q) 



-ej exp 



ap 

3nm 2 



(3Xo 



ln(/3A) 



lA{p,q) , \p\ < m. 



(39) 



(40) 



(41) 



As was already mentioned in the Introduction, the question as to whether the infrared divergence of the effective field 
matters at large distances can be resolved only after summing up the infinite series of local terms proportional to the 
(5-function derivatives. We now see that the soft-photon contribution leads to vanishing of the effective electric field at 
all distances from the charge: substituting the above expression for J(p, q) into Eq. ([5]), and taking the limit Ao — > 
yields 



Af(x)=0 for all x. 



V. DISCUSSION AND CONCLUSIONS 



We have investigated the effective electric potential of a localized electron using two different regularizations of the 
infrared divergences. In both cases, our results suggest that the effective field calculated using the perturbation theory 
vanishes. Though leading to the same conclusion, the two regularization schemes are quite different. While vanishing 
of the field in the A-regularization resembles that of the S'-matrix amplitudes with a fixed number of real photons, and 
occurs as the result of summation of all divergent diagrams in the infinite perturbation series, the field nullification 
in the e-regularization takes place already at the one-loop level. To the best of our knowledge, nothing similar has 
been encountered in calculations of other field-theoretic expectation values. In addition to that, this scheme reveals 
a number of important cancelations between various diagrams, described in Sec. IIIII 

It should be stressed, however, that our consideration is inherently perturbative, and hence is not complete, because 
the power-expansion with respect to the coupling constant breaks down, explicitly or implicitly, in both schemes. As 
to the e-regularization, we found that this scheme is quite natural in view of the fact that the finite-temperature effects 
provide an effective infrared regulator proportional to the photon bath temperature, which renders any diagram with 
a given number of loops infrared-finite. But the resulting mixing of orders in a makes it necessary to take into account 
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diagrams with arbitrary number of loops. This complication bears some resemblance to that caused by the gluon 
mass generation in high-temperature QCD [281 ] . the difference being that in our case, it is imaginary contribution to 
the mass that is responsible for the breakdown of the perturbative expansion. 3 

On the other hand, as we showed in Sec. HVi the A-regularization is free of this drawback. This scheme allowed us 
to consistently take the limit e —> in the cut integrals, and rendered diagrams with 2-vertices vanishing. Specifically, 
the assumption Ao 3> permitted us to overcome the question about generation of an effective e, which is 

formally zero in this regularization. In itself, this does not contradict the fact that it is nonzero in the other scheme, 
because, as was mentioned above, the perturbative calculation using the e-regularization is not complete. However, 
appearance of an imaginary part in the particle mass seems to be a general property of thermal field theor y, signify ing 
the presence of a relaxation mechanism driving the quantum field to the equilibrium with the heat bath [25l l26l l29( . 
Therefore, if the result E is scheme-independent, the limit Ao — > cannot be taken in Eq. Pdj) , because in effect, 
it would violate the assumption Ao 3> |s|/m. Rather, Eq. (|41[) can be used in this case to estimate the function J(p, q). 
Indeed, the only formal condition that the parameter Ao must satisfy is Ao <C T (see Sec. IIV Al . At the same time, 
E/m ~ aT <C T, so that a rough estimate for J(p, q) can be obtained by replacing Ao — > E/m. Equation (|4T|) then 
would give J ~ — e7° exp |O(a )} I\, showing that if E =/= 0, the consideration based on the A-regularization goes 
beyond the perturbation theory, too. 

We would like to emphasize that this conclusion about essentially non-perturbative nature of the infrared problem 
is directly related to the question of whether the infrared divergences affect the long-range behavior of the effective 
field. We have seen that despite the infrared-divergent terms in the effective field are local at every loop order, the 
sum of the infinite number of such terms is not (see Sec. IIV B]) . Specifically, if the limit Ao — > can actually be taken 
in Eq. (f^T|) . the effective field vanishes at all distances from the charge. On the other hand, if the temperature effects 
set a lower 0(a)-bound on Ao, then the infrared factor in J(p,q) corrects the effective field at large distances; though 
of the relative order 1/r 2 , this correction cannot be determined within the perturbation theory. 
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APPENDIX: WARD IDENTITIES AND A RELATION BETWEEN 2-TERMS AND E. 




FIG. 9: General structure of diagrams representing the function G^.^. The blobs denote irreducible parts - the photon 
polarization operator, II, and the proper electron-photon vertex, V. Double lines are the exact propagators. 

As we have seen, the perturbation expansion in powers of the coupling constant breaks down in the case E ^ 0. In 
view of this breakdown, it is natural to look for non-perturbative means for investigating the infrared problem. One of 
such means is the use of the Ward identities. Unfortunately, the introduction of e-regulator, necessary in the real-time 
techniques to deal with the products of singular functions, violates the usual Ward identities. Still, consideration of 
the unregularized identities makes sense by the following reason. As we have seen, the finite E plays the role of an 
effective e- regulator, though it necessitates analytical continuation to the negative values of e. Therefore, under the 
assumption that the exact electron self-energy admits this continuation, the existence of E ^ suggests that e can 
be consistently set equal to zero. 



3 That the thermal field theory in its present formulation is inherently non-perturbative was discussed from a more general standpoint in 
Ref. El- 
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To derive the Ward identities for our model, let us consider the following exact Green function 

G». k) (,x,y, z) = Tr (T c [^(z)J^(x)My)] e~^g) , (42) 

where J^(x) = ^(x^^tp^x) is the electromagnetic current, and the subscript indices indicate position of the corre- 
sponding time argument on the contour C (see Sec.|TI|). As before, the components relevant to the infrared properties 
of the effective field have j = k = 1 (other components can be dealt with using complex conjugation). Using the 
current conservation d^J^ix) = 0, we find for their 4-divergences 

^Tr (T c [M*)J?(*)Mv)] = ~ V)^ ( T c [M*)Mx)\ e~ m *Q) 

-5\x - z)Tr (T c [Mx)Mv)] e ~ P " 4 ' O) > ( 43 ) 



— Tr (T c [M*)J?(x)M)] e- f)H *Q) = S 4 (x - y)Tr (T c [^ 2 {z)Mx)} e^ H * e ) . (44) 

Diagrams representing GvV fe % are shown schematically in Fig. [9l To compute the left hand sides of Eqs. (|43|) . ([44]) . we 
shall use the Ward identity for the photon polarization operator IP" (a; — y), [20] 

d M n%(x-y) = 0. 

It follows from this identity that diagrams in Fig. GJa) do not contribute to the left hand sides of Eqs. (|4"5)l . (I44|) . 
Going over to the momentum space, these equations thus yield 

D< U > (q + P)p^ ln) (q + p, q)T>i^ (q) = (q + p) - {q) , (45 ) 



2)( 21 > + p)rf ul) (g + p, g) + D( 22 ' (g + p)I* n) (g + p, g)j !»(") (g) = 2)( 21 ) (g +p), (46) 



where and are the exact electron propagators and proper vertices. Notice that all other vertex components 

can be derived from these by complex conjugation. The first of these identities can be rewritten as 

p„rf m) (g + p,g) = (^(g))- 1 - (©^(g + p))" 1 . 

Using this in Eq. ([46]) . and taking into account that T>e\q) = 2) e (?) + I>i 22 ^(g) for go > 0, we bring the other 
Ward identity to the form 

p Ai r^ 211) (g+p,g)^(2)( 22 )(g+ P ))- 1 + (D( 11 )(<Z+p))- 1 . 

This identity is of particular importance as it gives insight into the connection between the appearance of E ^ and 
the vanishing of 2-terms, underlying factorization of the infrared contributions to the effective field. Namely, rewriting 
it as 

PM r (2ii) (l +P> 9) = 2iImS F (g + p), 

we see that r^ 211 j(g + p) cannot vanish, if E =/= (in fact, the exact vertex is singular for p — ► in this case). 
Accordingly, this function was found non-vanishing in the one-loop calculation of Sec. IIII CI On the other hand, it 
was demonstrated in Sec. IIV 51 that using the A-regularization, one can consistently set e = at any given loop order, 
and prove that the sum of diagrams involving 2-vertices vanishes. Accordingly, E is formally zero in the A-regularized 
perturbation theory. As discussed in Sec. [V] if E does not vanish in the complete theory, this result is valid only for 
A > E/m. 



[1] E. S. Fradkin, Trudy FIAN 29, 7 (1965); R. Fukuda and T. Kugo, Phys. Rev. D 13, 3469 (1976); R. Kallosh and 
I. V. Tyutin, Yad. Fiz. 17, 190 (1973) [Engl. Trans. Sov. J. Nucl. Phys. 17, 98 (1973)]; P. M. Lavrov, I. V. Tyutin, and 
B. L. Voronov, Yad. Fiz. 36, 498 (1982) [Engl. Transl. Sov. J. Nucl. Phys. 36, 292 (1982)]. 



Preliminary version - October 30, 2009 



17 



[2] N. Bohr and L. Rosenfeld, Kgl. Danske Vidensk. Selskab., Math.-Fys. Medd. 12, 3 (1933). 
[3] N. Bohr and L. Rosenfeld, Phys. Rev. 78, 794 (1950). 

[4] B. S. DeWitt, The quantization of geometry, in: Gravitation: an introduction to current research, ed.: L. Witten (Wiley, 
New- York, 1962) p. 266. 

[5] F. Bloch and A. Nordsieck, Phys. Rev. 37, 54 (1937); D.R. Yennie, S.C. Frautschi, and H. Suura, Ann. Phys. (NY), 13, 

379 (1961); S. Weinberg, Phys. Rev. 140, B515 (1965); G. Grammer and D. R. Yennie, Phys. Rev. D8, 4332 (1973). 
[6] T. Kinoshita, J. Math. Phys. 3, 650 (1962); T. D. Lee and M. Nauenberg, Phys. Rev. 133, B1549 (1964); G. Sterman and 

S. Weinberg, Phys. Rev. Lett. 39, 1416 (1977). 
[7] J. F. Donoghue, Phys. Rev. Lett. 72, 2996 (1994); Phys. Rev. D 50, 3874 (1994). 
[8] T. Toyoda, Phys. Rev. 77, 353 (1950). 
[9] Y. Nambu, Prog. Theor. Phys. 5, 614 (1950). 
[10] K. Hiida, M. Kikugawa, Prog. Theor. Phys. 46, 1610 (1971); K. Hiida, H. Okamura, Prog. Theor. Phys. 47, 1743 (1972); 

H. Okamura et al., Prog. Theor. Phys. 50, 2066 (1973); K. Yokoya et al, Phys. Rev. D 16, 2655 (1977). 
[11] R. Serber, Phys. Rev. 48, 49 (1935); A. E. Uehling, Phys. Rev. 48, 55 (1935). 

[12] I. Bialynicki-Birula and Z. Bialynicka-Birula, Quantum Electrodynamics (Oxford, Pergamon, 1975); J. Z. Kaminski, 

J. Phys. A: Math. Gen. 16, 2587 (1983). 
[13] E. P. Tryon, Phys. Rev. Lett. 32, 1139 (1975). 

[14] J. F. Donoghue, B. R. Holstein, and R. W. Robinett, Ann. Phys. 164, 233 (1985). 
[15] T. Altherr, Ph ys. Lett. 262B, 314 (19 91). 
[16] J. Manjavidze, |arXiv:hep-ph/951025i"j 

[17] H. A. Weldon, Phys. Rev. D 49, 1579 (199 4); Nucl. Phys. A566, 581c (1994). 
[18] D. Indum athi, Ann. Phys. 263, 3 10 (1998); |arXi"v:hep-ph/9607206| 
[19] A. Muller, |arXiv:hep-th/9912240| 

[20] N. P. Landsman and Oh. G. van Weert, Phys. Rep. 145, 141 (1987). 
[21] L. V. Keldysh, Sov. Phys. JETP 20, 1018 (1964). 
[22] J. Schwinger, J. Math. Phys. 2, 407 (1961). 

[23] A. J. Niemi and G. W. Semenoff, Ann. Phys. 152, 105 (1984); Nucl. Phys. B230 [FS10], 181 (1984). 
[24] N. P. Landsman, Nucl. Phys. A525, 397c (1991). 

[25] R. L. Kobes and G. W. Semenoff, Nucl. Phys. B260, 714 (1985); B272, 329 (1986). 
[26] W. Keil, Phys. Rev. D 40, 1176 (1989). 

[27] S. Weinberg, The Quantum Theory of Fields (Cambridge Univ. Press, 1995) Vol.1, Ch.13. 

[28] A. Linde, Phys. Lett. 93B, 327 (1980); D. J. Gross, R. D. Pisarski, and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 
[29] H. A. Weldon, Phys. Rev. D 28, 2007 (1983). 



Preliminary version - October 30, 2009 



